Synthesis of arbitrary, two-mode, high visibihty A^-photon interference patterns 
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Using coherent states and linear optics, we demonstrate the synthesis of arbitrary interference 
patterns and establish that neither the shape nor the visibility of A^'-photon interference patterns 
can be used as a quantum signature in general. Specific examples include saw-curve and rectangle 
curve interference patterns, and phase super-resolution with period shortening of up to 60 times 
compared to ordinary interference. The former two with a visibihty close to 100% and the latter 
with a visibihty in excess of 57 %. 
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The rapid development of experimental techniques has 
led to the demonstration of many remarkable quantum 
interference effects. A decade ago, it was shown that 
with specific A^-photon quantum states, one could break 
the Rayleigh diffraction limit and make optical interfer- 
ence patterns with a smallest feature size A^-times smaller 
than with ordinary light [1 . This discovery laid the foun- 
dation for quantum lithography, and as the name con- 
veys, it was thought that this was manifestly a quantum 
feature. However, superpositions and interference also 
manifest themselves in the classical world. It is therefore 
of interest to delineate what interference effects belong 
to the realm of the classical world, and which require 
quantum states. 

In different contexts limits arise for how large a visi- 
bility one can obtain using classical and quantum light. 
For instance, in the two-photon Hong-Ou-Mandel exper- 
iment, one can in principle achieve 100 % visibility with 
with both the input state 1 1) (g) 1 1) and the state 1 1) (g) |a) , 
where |a) is a weak coherent state. However, two mu- 
tually phase randomized classical input states will never 
reach a visibility in excess of 50 % [2]. Likewise, letting 
two classical states interfere in a Mach-Zehnder inter- 
ferometer and measuring the probability of detecting m 
and N — m photons respectively in the two output ports 
will likewise never result in a visibility of the A/A^-pcriod 
fringes > 50 % |3]. For three and four photon visibility 
experiments other limits to the obtainable visibility are 
81.8 % and 94.4 %, respectively [Hll]. However, all these 
limits are derived with the reference to a specific setup 
and a specific detection method. Here, we consider the 
interference between two coherent states using a general 
A^-photon projection measurement (realized with linear 
optics and coincidence detection). In this case we show 
that there is no difference between classical and quan- 
tum states neither regarding the visibility (that can be 
100 % in both cases) nor in the obtainable shape of the 
interference curve. Only by restricting the measurements 
to \m, N — m){m, N — m|-projectors one will regain the 
results in [3]. The distinguishing quantum signature is 
instead in success probability. Only with quantum states 



can one surpass limits to, e.g., phase sensitivity of classi- 
cal Hght [HIT]. This is regardless of the detection method. 
A general, two-mode, A^-photon state can be written 

N 

IV') = ^c„|n,Ar-n) 

r!=0 

N 

= 5]6„(at)"(6t)A'-"|0,0), (1) 

where 6„ = c„/ n\(N — n)!. If we formally divide 
X;fo«(a'^)"(&'^)^^" with (S^)^ and make the substitution 
{a} /w)" = z" we get the complex polynomial 

bnz" = bwiz - zi){z - Z2) • • • (z - zn), (2) 

Tl = 

where we have used the fact that any complex A^-th de- 
gree polynomial has exactly N complex roots. Hence, it 
is always possible to express any two-mode, A^-photon 
state 

1^) = 7j^^fiia^-^ib^)---^fN{a^~ZNp)\Q,0), (3) 

where A/'„, are real numbers such that Af^i^ + \zn\'^) = 1- 
Thus, any two-mode, A^-photon state can be written 
as a direct product of two-mode single photon states. 
Hofmann ^ made the important observation that to 
make a probabilistic projection measurement onto \4>) {ip\, 
one could split the state to be measured into A^ two- 
mode paths and using linear optics, unitarily transform 
the state in path n as A/'„(a'f — z„&1')|0,0) — )■ ajj|0,0). 
Then, in the case one detects one photon in the a„ 
mode of each path, one has in fact projected the in- 
put state onto \ip){'tp\. It is convenient to use polar- 
ization states and, e.g., take the two modes a„ and 6„ 
to be horizontally and vertically polarized modes. Then 
any single photon polarization transformation can be 
achieved with two wave plates: The first compensat- 
ing the relative phase between the a„ and 6„ mode by 
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6n = — Arg(z„) so that the (in general elhpticahy po- 
larized) state A/'„(aJj — z„63^)|0,0) is transformed to the 
linearly polarized state A/'„(a|( — |z„|&J^)|0, 0). The second 
being a polarizer rotated to the angle = arctan |z„| so 
that the linearly polarized state passes it, but the orthog- 
onal polarization is blocked. 

This method was further developed in [3] for measure- 
ment projection of NOON states, and demonstrated in 
[inilll] to project out the NOON state components from 
coherent states to demonstrate A/TV-period interference. 
In [12] the method was shown to provide a means of 
breaking the Heisenberg limit in intcrferometry when 
used with a \N/2, N/2) state input. A similar method 
was used in [131 HI! demonstrate resolution beyond 
the Rayleigh limit using independent light sources, such 
as thermal light. A scheme to synthesize arbitrary filter 
functions from N cascaded polarization "components," 
based on a similar factorization was suggested jl5j in 
1967. However, in this proposal the individual compo- 
nents were frequency modes and not single photon states. 

The stated factorization method can be taken further. 
Suppose the input state is a two-mode, linearly polarized 
coherent state 



oo oo 



|a,a) =exp(-|an ^ ^ ^ 
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(4) 



where for simplicity we take a to be real. We want 
to detect the interference between the phase-shifted in- 
put state and the projector where is given by 
Eq. ([T]). The unitary, differential phase-shift operator is 
U{(f)) = exp[-i(/)(a1'a - 6l'6)/2]. The detected TV-photon 
coincidence probability then becomes 

P(<^) (X \{a,a\U\m)? 



a 



a 



2N 



exp(-2a2) 



N 



exp(-2a^) 



^6„exp(i0[2n- 7V]/2) 

N 

K exp(i0n) 



(5) 



The sum within the absolute sign can be identified as 
a truncated Fourier series. Therefore, up to the high- 
est "frequency" iV0/2 in the parameter 0, any function 
can be expanded in the exponential function basis. This 
means that we can mimic any A^-photon state's projec- 
tion onto any A^-photon projector with a coherent state 
input! 

However, splitting a polarization, two-mode coher- 
ent state [/(0)|a, a) spatially into A^ paths results 
in a product state of A^ identical coherent states 
U{4>)\a/^/N,a/\/N). The fact that this is a product 
state makes each state's projection probability onto a 
certain single photon projector, v.i.z. 




FIG. 1: Projection measurement setup. A weakly excited 
two-mode coherent state, linearly polarized at 45 degrees from 
the vertical, impinges from the left a birefringent wedge im- 
parting the differential phase shift (j>. It is subsequently di- 
vided by a sequence of non-polarizing beam splitters into A*' 
paths. In each path, a certain single photon state is projected 
out by preceding the detector with a polarizer set at Qn and 
a birefringent wedge set to the differential phase On- Finally 
the detection coincidences are recorded. 



statistically independent of any other coherent state's 
projection probability [16] . The total probability is sim- 
ply the product P((^) = n^gP„((/)) of the A^ individual 
projection probabilities. Thus, if one uses coherent input 
states it is not necessary to measure the "clicks" in each 
path in coincidence. Instead, each single photon projec- 
tor n can be measured separately using a weak coherent 
state. The result Pn{4>) is recorded, and the final prob- 
ability P{4>) is subsequently obtained by multiplying all 
the individual probabilities as demonstrated for a NOON 
state in |11| . Hence, we need only to implement one of 
the arms in Fig. [T] at a time, and reconfigure the bifre- 
fringence 0„ and polarizer angle £>„ between each run. 
Note that this is only possible for coherent states. If, 
e.g., quantum states or thermal states are used, as is the 
case in [T^] and [TJ], respectively, there are correlations 
between the different single photon projection probabil- 
ities and a coincidence measurement is a must, whereas 
for coherent states either method works. 

To demonstrate the generality of the method, we first 
synthesized a rectangle function interference pattern. 
That is, we would like to have 



\{a,a\U{(j))\i;s)\'^ oc Rect((/),7r/2) = 



|0|<7r/2, 

1 otherwise. 



P„(0) = |(a/^,a/^|J7t(0)A/;(at_^Jt)|o,o)|2 (6) 



(7) 

To this end, for A^ = 30 which gives a reasonable 
rectangle-like function, we compute the 31 lowest Fourier 
expansion coefficients of the Rect function and identify 
the coefficients with the expansion coefficients 6„ in Eq. 
([5]). Subsequently the expansion is factored in single pho- 
ton projector terms by Eqs. ^ and (|3| (see Supplemen- 
tary Material). The thirty single photon projectors are 
then measured and multiplied. The result is displayed in 
Fig. [2] where the only fitting is the amplitude of the the- 
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FIG. 2: Rectangular interference pattern. The solid (black) 
line (215 data points) shows raw data and the dashed (red) 
line shows the theoretically expected curve with the amplitude 
chosen for best fit. At selected data points, error bars show 
the ±cr statistical uncertainty. 
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FIG. 3: Saw-curve interference pattern. The solid (black) 
line shows raw data and the dashed (red) line shows the the- 
oretically expected curve with the amplitude chosen for best 
fit. At selected data points, error bars show ±a statistical 
uncertainty. 



oretically predicted coincidence probability. The exper- 
imental curve represents raw data, no background sub- 
traction or any other data processing have been done. 
The maximum count rate employed in the experiments 
was around 5.2 MHz. A discussion about the error bars 
will follow below. 

Next, we synthesized a Saw-curve interference pattern. 
The interference pattern being the square of the overlap 
between the input state and the projector, we want the 
(27r-periodic) overlap between the input state and the 
measurement projector to be 

\{a,a\Um)\^\<l>M^/^ (8) 

in the interval {— 7r,7r}. Again finding the Fourier co- 
efficients for N = 30, plugging them into Eq. ([2]), and 
factorizing the polynomial to obtain the single photon 
projectors (see Supplementary Material), we finally get 
the raw data displayed in Fig. [3] where the only fitting is 
the amplitude of the theoretically predicted curve. This 
interference pattern has the particular feature that its 
derivative, which governs the phase sensitivity, is almost 
constant over large intervals. 
Finally, a NOON state 

IV'noon) = ^(|A^,0)-|0,A^», (9) 

can give an interference pattern with the smallest feature 
size ~ X/{2N), where A is the wavelength of the light [1]. 
This feature size reduction can be explained in terms of 
maximizing a state's dynamical evolution speed |17H20| 
or in terms of A^-photon quasi-particles having N times 
the linear momentum of the photons making up the quasi 
particles, and thus having a de Broglie wavelength N 
times smaller [2m24j . Up to now, there have been many 
proposals and demonstrations H [HHIIl HIMH] of NOON- 
state interference with = 2 to TV = 30. The A^ = 6 to 



A = 30 experiments |101 111) , were made with coherent 
light with the same method we now generalize. There has 
also been other demonstrations and proposals of multi- 
photon (i.e., using coincidence measurements) interfer- 
ence using non-entangled, classical states also showing a 
period shortening, but at the expense of a reduced inter- 
ference visibility. In, e.g., [3] periods down to 1/4 of the 
"regular" interference period were measured, but with a 
visibility of only about 29 %. Using the method above, 
we obtained interference patterns for projection of a co- 
herent state onto A^ = 10, 15, 30, and 60 NOON state 
projectors, shown in Fig. |4] We see that for A^ = 10 
we get an excellent fit with the predicted pattern. The 
poorest visibility, where the visibility of a fringe is de- 
fined by the adjacent local minimum and local maximum 
of the interference pattern is 99.5 % in this case. Even 
for A^ = 60 we get a reasonable fit with the theory and a 
minimum and maximum visibility of 57.5 % and 86.2 %. 

The visibility is limited by five effects. The first one 
is the smallest step with which we can vary (jj. It is set 
by the combination of birefringent material (quartz in 
our case), the birefringent wedge angle which is 19 deg, 
and the linear motorized stage minimum step size 1 /xm. 
This leads to a minimum phase difference resolution of 
29 mrad which for A = 60 is too poor (one period of 
the pattern, occupying 27r/60 radians, is probed by less 
than four measurement points). As seen in Fig|4](d), we 
simply "miss" many minima and maxima of the pattern, 
yielding a lower visibility than that set by other measure- 
ment errors and noise. For A^ < 30, the dominant source 
of error is the stochastically varying quantum efficiency 
of the detector. The measured coincidence probability is 
P(0) = n^^]^77„P„(0), where 7y„ is the quantum detector 
efficiency for the single photon projector n. In our case, 
we use the same detector for each projector, and therefore 
Vn where r = 480 s is the time it took between 
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FIG. 4: NOON-state coincidence patterns. In the figures on 
tlie left, tlie solid (black) lines represent raw data points con- 
nected by straight lines. No background subtraction has been 
done. The dashed (red) curves are the expected Asin'^ {N(f>/2) 
with A chosen for best fit. The figures on the right show a 
magnified part of the data points with error bars points and 
the fit curve. Error bars show ±a statistical uncertainty, a, 
iV = 10. h, N = 15. c, iV = 30. d, N = 60. The boxed area 
shows the portion with highest visibility, with the maximum 
of 86.2 % occurring at 1.75 tt. 



measuring the point (/) of one projector to the same point 
for the next projector. Assuming that the variation in 
quantum efficiency Ar](nT) = ri[nT) — (77(t)) is uncor- 
related between one projector to another, the expected 
standard deviation due to this error is 



((AP(<^))2)V2 = ViVF(^) ^M^^^)!^. 

V 



(10) 



This equation shows two important features. Firstly, that 
due to the multiplication of detection events the error 
will grow with the square root of the projector's photon 
number. This sets an upper limit to the method's appli- 
cability. Secondly, the absolute error is proportional to 
the coincidence probability, or stated otherwise, the rel- 
ative error is constant. This is clearly seen in our data as 
an excellent fit between experiment and theory as long as 
P{4>) is small. Another error source, but less important 
for our measurements, is small errors in the phases 0„. 



For a rapidly oscillating function, the exact position of 
each single photon projector's minimum has a substantial 
impact, as predicted in [29]. In addition, for the NOON 
state interference, one should ideally be able to get all 
the light in and out of one linear polarizer. At best, we 
manage a perpendicular polarization suppression ratio of 
46 dB, and typically the suppression rate is > 40 dB. 
Finally, the detector dark count also sets a limit. For 
our NOON state interference this error source is negligi- 
ble as the dark count rate is < 2 Hz and the maximum 
count rate is 1.02 MHz. The detector efficiency (if it were 
constant) does not matter in the experiment as the de- 
tected state is a coherent state irrespective of n and 0. 
A low detector efficiency (about 20 % for our detector) 
can easily be offset by a higher input intensity. 

With the Hofmann-Guo method [H [5], it is possible to 
synthesize any interference pattern to, in principle, the 
precision given by N. The method does not limit the 
obtainable visibility, and we have shown that in practice, 
even N — 60 NOON-state interference patterns can be 
recorded with > 50 % visibility. We have also demon- 
strated that unusual (non-sinusoidal) interference pat- 
terns can arise from multi-photon interference of coherent 
states. Hence, it is neither the shape nor the visibility of 
multi-photon interference patterns that delineate a gen- 
eral border between classical and quantum interference, 
although in special cases limits do apply. 

However, with a coherent state input it is not possible 
to get quantum phase sensitivity (i-e. break the 

standard quantum limit [30j [3T]). On the contrary, as 
discussed in some detail in [llj , for a given mean photon 
number, e.g., the NOON state interference curves have a 
lower phase sensitivity than ordinary (single fringe) in- 
terferometry with a coherent state. 

It is possible to directly scale up the demonstrated 
method to the classical regime using ordinary photo de- 
tectors, and perhaps it is here the method will find ap- 
plications. One reason we have stuck to working at the 
single photon detection level is to show that the coinci- 
dence method is very flexible in that allows the synthesis 
of any two-mode, multi-photon projector, albeit at the 
expense of exponentially decreasing probability of coin- 
cidence events with increasing N. Another reason is to 
explore what practical limits (i.e., how the influence of 
various sources of errors) one encounters when one scales 
single-photon coincidence measurements to large num- 
bers. 

When using coherent state input, the switch from ob- 
taining the interference curves by coincidence measure- 
ments to multiplication of single photon projection prob- 
abilities allows considerable savings in time and equip- 
ment. In our case, the maximum probability of de- 
tecting a photon in a single temporal mode, defined by 
the response time of the detector, is about 0.1 for each 
single photon projector. Hence, as the TV single pho- 
ton projectors give statistically independent results, the 
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probability of detecting, say, 10 photons in coincidence 
would be around 0.1^" = 10~^°. With our detectors 
with a response time of 45 ns it would thus take at least 
45 • 10~^ • 10^" = 450 s to get a single coincidence click, on 
average. Detecting each (temporal) mode sequentially, 
instead of detecting TV (spatial) modes in parallel, it took 
about 900 s to measure the 60 projectors for a single 
phase shift (/). Thus, taking 215 steps to cover the range 
of phase shifts from to 2tt radians, a whole = 60 in- 
terference curve was obtained in about 8 h which is long, 
but feasible. 

The authors thank C. Kothc, S. Takeuchi, 
K. Edamatsu, and S. Harris for helpful discussions. 
The work was supported by the Swedish Research 
Council (VR) through grant 621-2011-4575 and through 
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Supplementary material 
Methods 

We have implemented a multi-photon experiment us- 
ing a linearly polarized HeNe laser (whose polarization 
was further "cleaned" with a polarizer) with a power of 
15 mW, neutral density filters, two birefringent wedges 
mounted on step-motor drive translators (with a preci- 
sion of 1 /im) a polarizer mounted on a rotation stage 
(with a precision of 0.2 degrees), and a single photon sen- 
sitivity avalanche photo diode (APD), see Fig. [l] The 
intensity of the laser was adjusted by the neutral filters 
such that the mean time interval between two photons 
was 20 times longer (« 1/zs) than the dead time of the 
detector (45 ns). Thus, we ensure that essentially each 
recorded "click" stems from a single photon. The laser's 
polarization was carefully adjusted to be at 45 degrees 
from the horizontal, thus creating the desired input state 
I a, a), where a ^ 1 in the temporal mode basis of the 
detector. The calibration of the wedges was found to give 
a 27r phase shift with a relative translation of 0.215 mm. 

Following the method in [11 , each arm of the inter- 
ferometer was implemented sequentially and the results 
were recorded on the computer. For each arm n, the 
first birefringent wedge was set to impart the differen- 
tial phase shift <j> (from to 27r for each projector) on 
the input state, the second birefringent wedge imparted 
the differential phase shift 6'„ and the polarizer after the 
second wedge was rotated to the angle p„ to adjust the 
relative amplitude of the horizontal and vertical modes 
(see Fourier Expansion Coefficients). The state was sub- 
sequently focused on the APD and the number of "clicks" 
during one second was recorded as a function of in a 
computer memory. After all N single photon projector 
counts [proportional to P„ ((/)))] had been measured as a 
function of (j), the counts at each setting of (j) were multi- 
plied together yielding a function proportional to P{<j))- 

All the displayed data is raw, that is, without the 
subtraction of dark counts and it was obtained with a 
counter gate time of 1 s. 



Fourier expansion coefficients 

Consider a 27r-periodic rectangle function defined by 
the two conditions 



Rect((/>,7r/2) 



101 < 7r/2, 

1 otherwise. 



(11) 



and Rcct((/), 7r/2) = Rcct(0 + 2/7r, 7r/2) where / is an ar- 
bitrary integer. Since the Rect function is even, it is 
advantageous (and natural) to use an even number of 
projectors to synthesize the function implying that N 
should be chosen even. (The function can also be ex- 
panded for odd N but with a less pleasing result). The 



Fourier expansion of this function can hence formally be 
written 

N/2 

Rect(0, 7r/2) = ^ fe(2n+Ar)/2 exp(z(/)n) (12) 

n=~N/2 

where the expansion coefficients are 



5(2n+Af)/2 



1 

2tt 



Rect(0, 7r/2) exp{—i(f)n)d(j) 
1/2 71 = 0, 

_sm(nV2) ^ / Q 



(13) 
(14) 



Using Eq. ([I]), the associated polynomial for this state 
for, e.g., TV = 10 is hence 



z 

with the roots 



z'^ + z^ + z^ 1 + z 



10 



57r 



0, (15) 



Zl 


= Z*2 


= 0.967612-1- 


iO.252442, 


Z3 


— Zi 


= 0.723141 -1- 


i0.6907, 


Z5 


= 4 


= 0.313207 + 


iO.949685, 


Z7 




= -0.463687 


-t-iO.29332, 


Zg -- 


- 7* 

- zia 


= -1.54027-1 


- iO.974347. 



(16) 



The fact that the Rect function is an even function 
leads to the result that all the roots of the associated 
polynomial come in complex conjugate pairs. Insert- 
ing these roots into Eq. ^ it is evident how to im- 
plement the ten single projectors. For example, the 
n = 1 projector is implemented by introducing a bire- 
fringence of — Arg(zi) w —0.255 rad = —14.6 degrees be- 
tween the horizontal and vertical directions. This bire- 
fringence should be followed by a polarizer set at the 
angle arctan|zi| = 7r/4 rad. In Table |l] we list the set- 
tings of the birefringence and the polarizer angle for the 
10 projectors in degrees. It is seen in Fig. [2] that the 



Root n 


Qn 




1 


45.0 


-14.6 


2 


45.0 


14.6 


3 


45.0 


-43.7 


4 


45.0 


43.7 


5 


45.0 


-71.7 


6 


45.0 


71.7 


7 


28.7 


-147.7 


8 


28.7 


147.7 


9 


61.2 


-147.7 


10 


61.2 


147.7 



TABLE I: The experimental parameters for a TV 
function projectors. 



10 Rect- 



resulting function is not a perfect rectangle, but Gibbs 



7 



Root n 


Qn 


dn 


1 


56.2 


0.0 


2 


33.8 


0.0 


3 


66.6 


-74.0 


4 


66.6 


74.0 


5 


23.4 


-74.0 


6 


23.4 


74.0 


7 


21.7 


-141.7 


8 


21.7 


141.7 


9 


68.3 


-141.7 


10 


68.3 


141.7 



TABLE II: The experimental parameters for a = 10 Saw- 
function projectors. 



where FS is the Fresnel sine function. Numerically, the 
coefficients are: 





65 = 


2/3, 




= be = 


-0.1607, 


63 


= 67 = 


-0.0273, 


b2 


= bs = 


-0.0272, 


bi 


= 69 = 


-0.0109, 


bo-- 


= ^^10 = 


-0.0121. 



(19) 



The polynomial associated to the expansion is thus 



phenomenon makes the interference curve overshoot on 
the steep flanks. This effect can be reduced by the use of 
a Lanczos-Fourier expansion or a Cesaro approximation 
of the Fourier series, at the expense of getting less steep 
flanks in both cases. Going to higher photon numbers 
one could in principle reduce the wiggles in the interval 
{7r/2, 37r/2}, but the overshoot height would remain the 
same, only the width of the overshooting peak could be 
decreased. 

In a similar manner, an = 10 expansion of the 
Saw^/^ function 

N/2 

Saw^/^(0, 27r) = ^ 6(2„+jv)/2 exp[-z(/in]. (17) 

n=~N/2 

Computing the expansion coefficients in the same manner 
as before, one arrives at 



->{2n+N)/2 



/2nsin(n7r) - FS(V2n) 
\/27rn3/2 



(18) 



2z^ 



0.1607(z^ + z^) - 0.0273(z^ + z'^) 

-0.0272(^2 + z^) - O.Gll(z + z^) 

-0.012(1 + 2^") = (20) 



with the roots 



zi = 1.49215, 

Z2 = 0.670175, 

Z3 = yl = 0.637336 -I- i2.22359, 

Z5=y* = 0.119116 -HiO.415582, 

zr^y; = -0.311654 + iO. 245907, 

ZQ=ylo ^ -1.97752 + il. 56034. 



(21) 



The settings, in degrees for the birefringence and the po- 
larizer angle are given in Table [11} 



